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We report on the electronic properties of few interacting electrons confined in a parabolic quantum 
dot based on a theoretical approach developed to investigate the influence of Bychkov-Rashba spin- 
orbit (SO) interaction on such a system. We note that the spin-orbit coupling profoundly influences 
the energy spectrum of interacting electrons in a quantum dot. Here we present accurate results for 
the energy levels and optical-absorption spectra for parabolic quantum dots containing upto four 
interacting electrons, in the presence of spin-orbit coupling and under the influence of an externally 
applied, perpendicular magnetic field. We have described in detail about a very accurate numerical 
scheme to evaluate these quantities. We have evaluated the effects of SO coupling on the Fock- 
Darwin spectra for quantum dots made out of three different semiconductor systems, InAs, InSb, 
and GaAs. The SO coupling on the single-electron spectra manifests itself by primarily lifting of 
the degeneracy at zero magnetic field, rearrangement of some of the energy levels at small magnetic 
fields, and level repulsions at high fields. These are explained as due to mixing of different spinor 
states for increasing strength of SO coupling. As a consequence, the corresponding absorption 
spectra reveal anticrossing structures in the two main lines of the spectra. For the interacting 
many-electron systems we observed the appearence of discontinuities, anticrossings, and new modes 
that appear in conjunction with the two main absorption lines. These additional features arise 
entirely due to the SO coupling and are a consequence of level crossings and level repulsions in the 
energy spectra. An intricate interplay between the SO coupling and the Zeeman energies are shown 
to be responsible for these new features seen in the energy spectra. Optical absorption spectra for 
all three types of quantum dots studied here show a common feature: new modes appear, mostly 
near the upper main branch of the spectra around 2 tesla that become stronger with increasing SO 
coupling strength. Among the three types of systems considered here, optical signatures of the SO 
interaction is found to be the strongest in the absorption spectra of the GaAs quantum dot, but 
only at very large values of the SO coupling strength, and appears to be the weakest for the InSb 
quantum dot. Experimental observation of these new modes that appear solely due to the presence 
of SO coupling would provide a rare glimpse on the role of SO coupling in nanostructured quantum 
systems. 

PACS numbers: 71.70.Ej,72.25.Dc,72.25.-b 



I. INTRODUCTION 



Impressive developments in nanofabrication technol- 
ogy have made it now possible to design the quantum 
dots and coupled quantum dots at the nanoscale. These 
systems comprise of a few electrons that are quantum 
confined, for our present purpose, at the semiconductor 
interface to form the zero-dimensional systems. What is 
more remarkable is that the electronic states in these sys- 
tems can be precisely controlled via the external voltages 
0. Magneto-optical studies of parabolic quantum dots 
(described as artificial atoms by us in 1990) have been 
intensely explored for more than a decade 0, 0, S @ 
in order to understand its unique electronic and opti- 
cal properties and as promising candidates for optoelec- 
tronic devices, applications in optical quantum informa- 
tion technology, etc. 0>I3- As a result of these studies, 
a very good theoretical and experimental understanding 
of the single-electron states in the dot has already been 
achieved. At the most basic level, the solution of the 
Schrodinger equation for an electron confined by a har- 



monic potential, v c — ^m*ujor 2 , ujq is the confinement 
potential strength, in the presence of an external perpen- 
dicular magnetic field is well known since the beginning 
of quantum mechanics 0,0. The eigenvalues in this case 
are given by 

E nl = (2n + \l\ + l)Ml- \lhw c 

where n ~ 0,1,2,... and / = 0, ±1..., are the princi- 
pal and azimuthal quantum numbers respectively, fl 2 = 
[ujq + j\u>c] i an d w c is the cyclotron frequency. Dipolc- 
allowed transitions among these energy levels will have 
energies USQ 

AE ± =m± \huj c . 

This relation has been verified to great accuracy by a 
variety of experiments 3,0, El- Interestingly, however, 
the observed magnetic field dependent FIR absorption 
in quantum dots containing more than one electron was 
found to be essentially independent of the number of elec- 
trons confined and instead was dominated by the above 



2 



relation for AE± p|. It was a rather puzzling result be- 
cause according to this, magneto-optics was clearly inca- 
pable of providing any relevant information about the ef- 
fect of mutual interactions of the confined electrons. The 
puzzle was later resolved by Maksym and Chakraborty 
[H d j wri o pointed out that for a parabolic QD 
in an external magnetic field, the dipole interaction is 
a function of the center-of-mass (CM) coordinate alone 
and the inter-electron interaction does not play any role. 
Despite this somewhat disappointing performance of a 
parabolic dot, FIR spectroscopy of QDs (parabolic or 
otherwise) has generated enormous interest for over a 
decade that is yet to subside |6| . In this paper we include 
another interesting element into the problem, the spin- 
orbit (SO) interaction. We demonstrate here that in the 
presence of spin-orbit coupling the magneto-optical tran- 
sitions show many new and interesting features that can 
be tuned by the SO coupling. Based on these results, we 
propose that magneto-optical transitions are best suited 
to determine optically the unique effect of SO coupling 
in quantum dots described below 

Interest on the role of the spin-orbit coupling in nanos- 
tructured systems is now at its peak, due largely to its 
relevance to spin transport in low-dimensional electron 
channels [T2I Il3j . The intriguing possibility of tuning 
the SO field and thereby coherently manipulate electron 
spins in quantum dots has sparked major activities in 
the pastyears Ejlllll [13, E (13, Hfl [13, ISE E3, 
EllMliillllilWlIJ. It is hoped that an improved 
understanding of spin dynamics in the QDs might pave 
the way for future electronic and information processing, 
especially in quantum computing and quantum commu- 
nication [32|. Spin degree of freedom is perhaps more 
advantageous than charge because unlike charge, spin is 
not coupled to electromagnetic noise and therefore have 
much longer coherence time 33] . Improved knowledge of 
the influence of spin-orbit coupling in quantum dots is 
therefore quite essential in this pursuit. While the ma- 
jority of experimental efforts has focused on magneto- 
transport measurements |34j. here we present the results 
for the optical absorption spectra that are experimentally 
observable and could, in principle, provide an important 
probe of SO coupling in few-electron quantum dots. 

The spin-orbit interaction in semiconductor het- 
erostructures can be caused by an electric field perpendic- 
ular to the two-dimensional electron gas (2DEG). Riding 
on an electron, this electric field will be felt as an effective 
magnetic field lying in the plane of the 2DEG, perpen- 
dicular to the wave vector k of the electron. The effective 
Zeeman interaction of the electron spin with the field lifts 
the spin degeneracy (internal Zeeman effect). This is usu- 
ally referred to in the literature as the Bychkov-Rashba 
mechanism. This results in an isotropic spin splitting 
energy Aso at B — proportional to k ^ • 

Let us consider an electron in the 2DEG moving with a 
velocity v in the presence of an electric field E. In the rest 
frame of the electron, this transforms (relativistically) 



into an effective magnetic field -B c ff, 

^ cff = ~2c^ X ^ 

where c is the speed of light. The magnetic moment of 
the electron will then couple to B c ff . The resulting spin- 
orbit interaction is, 



1 



eft 2 



; a-(kx E) 



pj ■ [v x E) 



since, /i = a, and v = hk/rriQ. It can be rewritten 

2too 

as 



T~tso 



eh 2 



AitIqC 2 



a-{kxE)=a'{E z )(-iV xa) z , (1) 



where the electric field is aligned along the z axis. 

Alternatively, a general spin-orbit Hamiltonian that 
stems directly from the quadratic in v/c expansion of 
the Dirac equation is [3(| 



(2m c) : 



r W(r) ■ (a x <p). 



(2) 



The electric field associated with V{r) is E(r) — W(r), 
and is directed along the z direction. The spin-orbit in- 
teraction Hamiltonian 

eh 2 1 1 

^so = -77, ^)r(?xp)-n = a'(£ 2 )-(ffxp)-n 

{2mQc) z ft ft 

— a'(E z )(a x k) ■ h, 



where a 1 = e(H/2moc) 2 is then identical to that in 
Eq. An important point to note here is that a non- 
vanishing gradient in Eq. J3J) requires that the system 
must have inversion asymmetry. In the present case that 
arises from the structural inversion asymmetry |3?| . The 
spin-orbit interaction that we are here concerned with is 
therefore described by the Hamiltonian 



n 



so 



a(k x a) z 



la a 



d_ 

' dx 



d_ 

' dy 



(3) 



where the z axis is chosen perpendicular to the 2DEG 
(in the ccy-plane), a is the spin-orbit coupling constant, 
which is sample dependent and is proportional to the 
interface electric field that confines the electrons in the 
x ~ y plane, a — {a x , a yi a z ) denotes the Pauli matrices, 
and k is the planar wave vector. This is the Bychkov- 
Rashba Hamiltonian |3^ that has been receiving of late 
rather widespread attention [34|. 
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The single-electron Hamiltonian for the 2DEG includ- 
ing the Bychkov-Rashba term has the form 



n 



— h — (a x p) 

2m* h 1 yJz 

V 2 ia (a — -a — 

2m* \ v dx x dy 



-V 2 aV" 

h 2 - 



2m 



2m* 



where V 2 = d 2 /dx 2 + d 2 /dy 2 and V ± = d/dx± id/dy. 

Since the operators p x and p y commute with the 
Hamiltonian, we can search for a ^ eigenstates of the 
form 



k y ) = e ik * x+ik v y C °\ G ) = e lk * x+iky 



C+ 

c- 



with | cr) = ( s P m U P) or ^ x ) ( S P' H down). Solu- 

tions of 



are readily obtained as 



1 

71 



±fc„ =F ife I e 



ik x x-\-ikyy 



The energy dispersion then consists of two branches 

£±(fc) = —k 2 ±ak 
y ' 2m* 

with an energy separation Ago = £ + — £~ = 2afc for a 
given k. The spin parts of the wave functions X^{k x , k v ) 
are mutually orthogonal and (x ± l^z lx ± ) = 0- Therefore 
in the states "J/ the spins of the electrons lie in the xy- 
plane and point in opposite directions. In addition, 



(x ± l^lx ± > 



k 



(x ± l^lx d 



1k x 



i.e., the spins are perpendicular to the momentum 
(k x ,k y ). Spatial alignment of spins therefore depends 
on the wave vector |l2l l35l | . The Fermi surface is a pair 
of concentric circles with radii kp^ max and kp^min- In the 
present paper, we are dealing with systems having rota- 
tional symmetry. The formalism in that case is derived 
in detail in Sect. II. 

Several experimental groups 34] investigating the 
Shubnikov-de Haas (SdH) oscillations in a 2DEG con- 
fined at the heterojunctions with a narrow-gap quantum 
well (e.g., InGaAs/InAlAs, InAs/GaSb, etc.) have al- 
ready established that lifting of spin degeneracy results 
from inversion asymmetry of the structure which invokes 



an electric field perpendicular to the layer. Experimen- 
tally observed values of the SO coupling strength a lie in 
the range of 5 - 45 meV.nm [34|. Energy levels of two 
interacting electrons confined in a parabolic quantum dot 
in an external magnetic field were recently reported by us 
for this range of SO coupling strength [^H . In the absence 
of SO coupling, electron-electron interaction causes the 
ground state energy to jump from one angular momen- 
tum value to another as the magnetic field is increased 
0,0. The influence of the SO coupling is primarily to 
move the energy level crossings to weaker fields [3l|. 

Our paper is organized as follows: In Sect. II, we 
present the essential formalism for our study of non- 
interacting electrons in parabolic QDs in the presence 
of spin-orbit interaction. The single-electron basis and 
the dipole matrix elements are derived here. We also 
explain why the dipole-allowed transitions are so signif- 
icantly influenced by the presence or absence of the SO 
interaction. The classic Fock-Darwin spectra for three 
different quantum dots systems (In As, InSb, and GaAs) 
with or without the SO coupling are presented and dis- 
cussed in detail. The formalism for the many-electron 
system in a parabolic QD with SO coupling is presented 
in Sect. III. The complexities of introducing the SO cou- 
pling, in particular for interacting electrons, are made 
clear in Sect. III. The task of finding a suitable numeri- 
cal technique is even more challenging, and an approach 
that is appropriate for our purpose is described in the 
Appendix. Numerical results for the energy levels and 
the optical absorption spectra for the three types of QDs 
containing upto four interacting electrons are presented 
and discussed in Sect. IV. It should be pointed out that 
the low-lying energy levels calculated here for the single- 
and multi-electron quantum dots can, in principle, be ob- 
served in transport |3^ , or capacitance spectroscopy [3^| . 
Given the accute interest on the influence of SO coupling 
in nanostructured systems and the resulting intense ac- 
tivities on this topic, it is no surprise that many differ- 
ent theoretical techniques have been put forward in the 
literature. To view our work in proper perspective, we 
present a brief review of many of those theoretical papers 
in Sect. V. We conclude with a brief outlook for future 
work along this direction in Sect. VI. For a brief account 
of our earlier work on SO coupling effects in parabolic 
QDs, see Refs. [HIHI. 



II. SINGLE-ELECTRON PICTURE 

The Hamiltonian for an electron in a parabolic con- 
finement and under external magnetic field is given by 



Ho 



2m 
a 

n 



e -t)+\m^ 2 r 2 



G X \p 



e 



-.4 



1 



-gfi B Ba z 



(4) 
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Here a is the vector of Pauli matrices, i.e. 



a = a x i + <7 y j + a z k 
1 \ f , ( 



1 / l+ \ i 



We work in the symmetric gauge and the vector poten- 
tial corresponding to the external perpendicular mag- 
netic field is 



A = -(-y,x,0). 



(6) 



The term § 



a x 



in the Hamiltonian is the 



spin-orbit (SO) coupling due to the inhomogenous po- 
tential confining the electrons to the 2D plane and possi- 
ble external gate voltages applied on the top of the dot. 
The parameter a determines the strength of this coupling 
and, in case of external gate voltages its magnitude can 
be varied. Finally, the last term \g\iBBa z is the ordinary 
Zeeman coupling g being the effective Lande' ^-factor. 
The eigenstates of the single particle problem 



Hq4> = £0 
are clearly two-component spinors 

|A) = 



(J) 



(8) 



Writing the equation Q in polar coordinates and substi- 
tuting a trial wave function of the form 



/T(r)e 



we 



(9) 



it is easy to see that the quantum numbers ft and ft 
must be integers and that they depend on each other in 
the way 



ft = ft - i. 



(10) 



Hence we need only one quantum number for the angular 
motion, i.e. solutions of the single particle equation (JTJ 
are of the form 



|A) = |M> 



(11) 



Here the quantum number k is associated with radial 
motion (and not to be confused with the wave vector 
described in Sect. I). The form of the spinor (|ll|l sim- 
ply restates the fact that under SO coupling the good 
quantum numbers are related to L + S. In our case the 
conserved quantity is 



j = ft' 1 + si' 1 =t+\ 



(12) 



where s z = ±^ depending on the component of the 
spinor, i.e., + i for the upper component and — i for the 
lower one. 



In order to find the radial wavefunctions /I'l we trans- 
form to dimcnsionless units by setting 





eB 


LO c = 


* ' 




m*c 




r 2 


X = 










ea 2 B 


b R = 


h 2 ' 


hi __ 




i ~ 


4w (1 H 



m*uJo(\+uj 2 c lAujly ' 
m*aa 



h 2 ' 
m ± = i± 2b R , 



gfi B B 



2 — ±- 

Substituting these into the Hamiltonian (0J the radial 
part of the equation Q takes the form 



J" _i_ „T' 



ft x 



■>■!! - <! \ " ~ ^ - i + 4 J 3 T 

^ 1/2 (/ + ^/ + &«/)=0 (13) 
(l+l) 2 x t \ I 



fix 1 



1/2 „V _ _L„T _ h _j 



2x~ 



9 <-b R g<\=0 (14) 



of two coupled differential equations. For (3 = (i.e., 
a = 0) these equations describe radial motions of two 
independent two-dimensional harmonic oscillators with 
eigenvalues 



T.l 

vl't = n - 



+ 1 



(15) 



and the eigenfunctions 



gni 



(n + \e\y. 



-*/* x W* L W( x ). (16) 



Here L' n is the associated Laguerre polynomial defined, 
for example, by the formula 



LW(x)=e*x-W^-(e-*x n +W). 
n y ' dx n V J 

Therefore it is logical to seek the solution for Eqs. (|13l 
114(1 in the form of the expansion 



(17) 



n=0 



In spinor language this corresponds to the expan- 
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sion 



|A) = 



JIB 







n=0 

ol (£+i)e 



e c i,f+i 



n=0 



(x) 



(18) 



The coefficents can be obtained by minimizing the 
expectation value (A|Wo|A). At this point it is usefull to 
relabel our spinors. For example, for non-negative values 
of £ we set 



u 2n = e 



9n,t( X ) 





Z2n 
Z2n+1 





g n ,e+i(x) 



■-n,e 



1 



(19) 

(20) 

(21) 
(22) 



The minimization of 

<A|Wo|A> = £ 

n,n' 

leads now to the diagonalization of the symmetric tridi- 
agonal matrix TL with diagonal 

diagW = (i/Jj, . . . , i/^, ^ +1> . . .) (23) 

and subdiagonal 

subdiagH = — (0, \/£ + l?y + , 77", . . . , \fnq~ , 

Vn + £ + lr) + , Vn + It)-, . . .V (24) 

The diagonalization yields a set of solutions: the set 
f eigenvalues and the set {zn } of eigenvectors 
indexed by the particular solution k and the fixed angular 
momentum £. These are the energies of the spinors and 
the expansion coefficients (12 11221 . For negative values of 
£ the tridiagonal matrix consists of the diagonal 



1 T 



■) 



(25) 



diagW = {vo,t+i^l, 
and the subdiagonal 

subdiagW= -^(0, VW\v~, V + , ■ ■ ■ , Vnrj + , 

v/n+ |^|?7^,V^TT?7 + ,...). (26) 

The spinors thus obtained comprise our single-particle 
basis 

B s = {|Ai)|i = 0,l,2,...} 

= {|fc,*)|fc = 0,l,2,...;* = 0,±l,±2,...}. 

(27) 



A. Dipole matrix elements 

According to the Fermi golden rule the intensity of 
absorption in dipole approximation is proportional to the 
square of the matrix element 



I=</l£ne ±w «|i> 



i=l 

when the transition goes from the initial A^-particle state 
\i) to the final state |/). To evaluate this we need to 
know the dipole matrix elements d\',\ between the spinor 
states I A') and |A), i.e. the elements 

dv,A = (A'|re ±l9 |A) - {k\l'\r e ±lB %l) . (28) 

For simplicity we consider only the circular polariza- 
tion e +i0 (the other circular polarization is obatined by 
reversing the roles of A and A'). Substituting expansions 
(fP7|) into the above expression we get 



dx> 



A', A 



+ c'ic l n Vn + £ + 2- c'i^c^Vn (29) 
when £ > 0, 

d\>,\ = a^/,oE [ c 'l c lV n + 1 ~ c 'i+i c W n + 1 

n 

+ c 't c W n + 1 - c 'n-i c lVn (30) 
when £ = — 1 and 

dv,A = a8e',e+i E [ c 'l c W n ~ ~? ~ c '1+i c W n + 1 



+ C 



/l n c l n Vn-£- 1 - c'^ +1 4'V / n + 1 



(31) 



when £ < — 1. The intensity is then obtained from T oc 
MaiA 2 1 2 E3- I 11 an our figures for the absorption spectra, 
the size of the points is proportional to the calculated 
intensity. 

The reason why dipole-allowed transitions in a 
parabolically confined quantum dot can be very different 
in the presence of SO interaction is explained as follows. 
When subjected to the radiation field with amplitude a 
and polarization e, the vector potential A in the single 
particle Hamiltonian 



Ho 



2m 
a 

+ H 



1 f-. e A 2 1 



*, ,2^.2 



m oj n r 



c 

a x I p 



-A 

c 



1 



g/.i B B(J z 



must be replaced with the potential 



A — > A + A U ,A U = eae 



ik-r—itut 



FIG. 1: Some of the low-lying energy states for non- 
interacting electrons confined in a InAs quantum dot for var- 
ious values the SO coupling strength (in meV.nm), a = [a], 
a = 20 [b], a = 30 [c], and a = 40 [d]. The labels in (a) are 
explained in the text. 



FIG. 2: Some of the low-lying energy states for non- 
interacting electrons confined in a InSb quantum dot for var- 
ious values the SO coupling strength (in meV.nm), a = [a], 
a = 20 [b], a = 30 [c], and a = 40 [d]. 



In the dipole approximation we assume that 
and correspondingly the Hamiltonian will be |4l| 
where 

In a many-body system when a = the first term gen- 
erates the CM density excitations where mutual inter- 
actions play no role. Consequently (in dipole approxi- 
mation) only transitions between these modes are pos- 
sible. When a is different from zero, the second term 
(oc a x e y — cr y e x ) in H! can create spin-density oscillations 
and interactions have effect on their properties. It is to 
be noted that, in SO coupled systems the dipole operator 

eCL 

still retains its familiar form, Q = — e-r, as is easily ver- 

c 

ified by evaluating its commutator with the Hamiltonian 
Ho 

[Q,H ] = iKH'. 

Dipole operator is independent of the electron spin. The 
dipole-allowed optical transitions are always between the 
same spin states, but the angular momenta must differ by 
unity. In the presence of SO coupling, neither the dipole 
operator nor the selection rule changes, but the SO inter- 
action mixes the neighboring angular momentum values 
(I and I + 1) as well as the spin and hence the selection 



rule now applies to the total angular momentum J as 
well. Therefore, transitions from other states that are 
not allowed without the SO coupling, are now allowed. 



B. Fock-Darwin spectra 

In our numerical investigations, we choose InAs, InSb 
and GaAs quantum dots with parameters, to* /mo = 
0.042, e = 14.6, g = -14, to*/to = 0.014, e = 17.88, g = 
-40, and to*/to = 0.063, e = 12.9,. g = -0.44, respec- 
tively. While the InAs quantum structures have been 
the system of choice for investigation of spin-related phe- 
nomena |34| , InSb quantum dots are interesting for their 
very high g values and a relatively large a (~ 14 meV nm) 
42]. For the GaAs quantum dots, the observed value of 
a is ~ 6 meV.nm [28|. In all these systems we consider 
the confinement potential strength to be fiioo = 7.5 meV. 
Some of the low-lying states of the Fock-Darwin spectra 
of the InAs, InSb and GaAs QDs are shown in Figs.Q]- 
13 respectively and the corresponding optical absorption 
spectra in these systems are presented in Figs.0]-|H| 

As compared to the Fock-Darwin spectra of quantum 
dots without the SO coupling (shown in panels Figs-^- 
El[a]) the most outstanding features in the energy spec- 
tra of quantum dots with SO coupling are the lifting of 
degeneracy at zero magnetic field, the rearrangement of 
some of the levels at small fields and level repulsions at 
higher magnetic fields. To get some insight into the mech- 
anism causing this kind of behavior, let us have a closer 
look, as an example, at the energy levels involved in the 
lowest absorption lines of the InAs dot [curves labelled 
"0" - "3" in Fig. [T] (a)]. 
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FIG. 3: Some of the low-lying energy states for non- 
interacting electrons confined in a GaAs quantum dot for var- 
ious values the SO coupling strength (in meV.nm), a = [a], 
a = 20 [b], a = 30 [c], and a = 40 [d]. 
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TABLE I: Schematic spinors corresponding to four Fock- 
Darwin levels, marked - 3 in Fig. (a) of an In As dot 
without SO coupling. The black rectangles stand for non- 
zero radial wave functions. 



In the absence of the SO coupling energies of these 
levels are given by the formula l|15(l . The corresponding 
spinors (schematic) are depicted in Table [I] where the 
numbers in the first column refer to the labels in Fig. 
(a). The spinors for electrons with negative and posi- 
tive Lande' ^-factors are shown in the middle and third 
columns respectively. In actual physical systems, conven- 
tionally only the spinors with g < are of any interest. 

The spinor states of electrons on lines and 1 of Fig.Q] 
(a) differ only by the orientation of the spin: on line 
the spin is parallel to the magnetic field (||z) while on 
line 1 the spin is antiparallel to the field. Thus the energy 
difference between these states is the Zeemen splitting. 
The total single-particle angular momenta (|I2[) are cor- 




FIG. 4: Optical absorption spectra for non-interacting elec- 
trons confined in a InAs quantum dot for various values the 
SO coupling strength (in meV.nm), a = [a], a = 20 [b], 
a = 30 [c], and a = 40 [d]. 



respondingly j — Since under the SO coupling j is 
a good quantum number these two states will never mix 
even when the SO coupling is on. When the coupling 
strength a increases, the higher lying states with 3 = \ 
couple to the state as well as states with j ' = — i couple 
with the state I. In Fig^this shows up as an increasing 
splitting of the lines and 1 when going through the pan- 
els from [a] to [d]. It should be noted, however, that the 
mixing has very minor effect on the ground state since 
the other states with j = \ are energetically very far 
from that. 

Turning now our attention to electrons on lines 1 and 

2 we see from the Table H] (column g < 0) that their 
spinor states both have the same angular momentum 
j = i + i = — |. Consequently the SO interaction can 
mix these states. The mixing is particularly pronounced 
when the states are nearly degenerate, i.e. in the vicin- 
ity of the crossing point of lines I and 2. At moderate 
coupling strengths this mixing leads to level repulsions 
as shown in panels [b] and [c]. When the coupling is 
very strong, energetically higher states with 3 = —\ also 
become important in the mixing, leading to the imper- 
ceptibility of the level repulsion in panel [d]. 

We also mentioned lifting of degeneracies and rear- 
rangements of energy levels as features of the SO cou- 
pling under small magnetic fields. Since we are inter- 
ested in the absorption the most important states for us 
are the ones which can be reached from the ground state 
[j = i) respecting the dipole transition selection rule 
Aj = ±1. These are the lowest states with 3 = —\ and 

3 = | corresponding to the lines 2 and 3 in Fig.^(a) and 
to the spionrs 2 and 3 (with g < 0) in Table HJ respec- 
tively, for vanishing SO coupling. Since as the SO inter- 
action is stronger the larger is the (angular) momentum 
of the electron, the spinor 3 possessing maximum orbital 
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FIG. 5: Optical absorption spectra for non-interacting elec- 
trons confined in a InSb quantum dot for various values the 
SO coupling strength (in meV.nm), a = [a], a — 20 [b], 
a = 30 [c], and a = 40 [d]. 
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FIG. 6: Optical absorption spectra for non-interacting elec- 
trons confined in a GaAs quantum dot for various values the 
SO coupling strength (in meV.nm), a = [a], a — 20 [b], 
a = 30 [c], and a = 40 [d]. 



angular momentum 2 is affected more than the spinor 
2. Hence, although the energies of both spinors are de- 
creased by the SO coupling the effect on the spinor 3 is 
larger. 

Keeping the above discusssions in mind it is now easy 
to interpret the features introduced by the SO coupling 
into the absorption spectra [Figs. JI]-|EJ]. Firstly, al- 
though the lower absorption branch consists mainly of 
transitions from the state to the state 2 it shows an 
anticrossing at moderate coupling strengths. This is a 
direct consequence of the mixing of the spinor states 1 
and 2 which results in two spinors, both with nonzero 
upper component. Thus we can see transitions from the 
ground state to both of these. Furthermore, as men- 
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FIG. 7: Some of the low-lying energy states for two interact- 
ing electrons confined in a InAs quantum dot and for various 
values of the SO coupling strength (in meV.nm) , a [meV nm] : 
(a) a = 0, (b) a = 20, (c) a = 30, and a = 40. 



tioned earlier the level repulsion between states 1 and 2 
resumes the form of level crossing when the SO interac- 
tion becomes very strong. This causes the anticrossing 
in the absorption spectra to disappear. In the InAs dot 
this happens already at a = 40 while in the InSb dot 
the anticrossing still persists. Since the Lande' ^-factor 
of GaAs is very small the Zeeman split state 1 does not 
meet the state 2 within the range of the magnetic field 
under consideration. Consequently we see at the lower 
right corners of Fig. HJ1 (b) - (d) only the beginnings of 
the lower branches of these anticrossings. 

Secondly, the upper absorption branch corresponds 
mainly to transitions from the state to the state 3 mod- 
ified by the SO coupling. The small magnetic field how- 
ever makes an exception. As we discussed above, at small 
fields the state 3 is energetically lower than the state 2 
due to the SO interaction. Thus we get a crossing of 
spectra at small fields. 

From Eq. (|15fl . the separation between the states 2 
and 3, and hence also the gap between the absorption 
line branches is roughly proportional to the cyclotron fre- 
quency lo c which in turn is linearly proportional to the 
magnetic field and inversely proportional to the effective 
mass of the electron. Thus, due to the very small effec- 
tive mass of the electron in an InSb dot the energy of 
the spinor 3 exceeds the energy of the spinor 2 already at 
very small magnetic fields. Consequently also the cross- 
ing of absorption lines of an InSb dot occurs at very small 
magnetic fields as can be seen in Figs. El(b)-(d). 

At this point it may be worth mentioning that the 
energetics of the single-electron quantum dot under the 
influence of the SO coupling and subjected to an exter- 
nal magnetic field depends strongly on the sign of the 
Lande' g-factor. This is contrary to the case without SO 
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coupling where the energy spectrum is independent on 
the sign of g although, of course, the orientation of spin 
is determined by it. Let us consider, for example, the 
lowest absorption branch in the case g > 0. From Ta- 
ble [I] we can deduce that also in this case the transitions 
mainly take the spinor to the spinor 2. Now, however, 
there is no spinor level which would cross or even come 
close to the energy of the state 2 and mix with it. Con- 
sequently the anticrossing described above would not be 
observable in this case. 

This concludes our discussion of the energy levels and 
optical absorption spectra in a non-interacting QD in a 
magnetic field and in the presence of SO interaction. In 
what follows, we describe the theory for a interacting 
few-electron parabolic quantum dot. 



III. MANY-ELECTRON SYSTEMS 

The basis Bn for N interacting electrons in a QD 
is constructed as a direct antisymmetrized product of 
single-particle basis Bs <f2*7)l of the form 

JV 

B N = ^0B s = {|A,)|i = l,2,...} 
i=i 

= A i2 ; . . .;X iN )\ij = 0,1,2, . ..} 

= {|fci,£i; . . . ; k Nl £ N ) \ kj = 0,1,...; 

4 = 0,±1,...}, (32) 

where A stands for the antisymmetrization opera- 
tor. It is also understood that the notations such as 
|Aj x ; Xi 2 ; . . . ; \ N ) represent the antisymmetrized direct 




FIG. 8: Dipole-allowed transition energies for two interacting 
electrons confined in a InAs quantum dot and for various val- 
ues of the SO coupling strength (in meV.nm), a [meV nm]: 
(a) a = 0, (b) a = 20, (c) a = 30, and a = 40. The size 
of the points in the figures is proportional to the calculated 
intensity. 
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FIG. 9: Same as in but for a three-electron InAs quantum 
dot. 



products, i.e 
|A, 



— I , Aj 2 , . 



; X tN ) =A[\X il ) <8 \X h ) 



■®\X h )]. 
(33) 

Usually it is possible to restrict the size of Bn using the 
conservation laws. For example, in a rotationally invari- 
ant system the total angular momentum is a good quan- 
tum number. Therefore, for example, we fix it to J, and 
accept into the basis only those states that satisfy 



N 



i=l 



The states of the interacting system are expressed as 
superposition of the non-interacting states taken from the 
basis set (|32|l 



i*> = E c ^ A ' i 



(34) 



To extract the coefficients Cj, we again resort to the min- 
imization, i.e. we minimize the Rayleigh quotient 



where TL is the total many-body Hamiltonian. Again 
this leads to the diagonalization of the Hamiltonian ma- 
trix with elements (Ai\H\Aj). The eigenvectors are the 
desired expansion coefficients and the eigenvalues the cor- 
responding energies of the interacting system. It will be 
clear from the Appendix that both these tasks, construc- 
tion of the Hamiltonian matrix and its diagonalization 
are numerically quite challenging. 
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A. Coulomb matrix elements 

We write the total Hamiltonian Ti as a sum of the 
single-particle operators J3J and two-body operators 
V(f,r f ) as 



N - N 



(35) 



Since our basis states | Ai) are diagonal by construction in 
Ho we only need to evaluate the matrix elements of the 
latter sum in (|33|) . Our many-body states | A^) are ex- 
pressed in occupation representation language (|33|1 . and 
therefore it is natural to proceed in the occupation num- 
ber space. This means that for the interaction part we 
have to evaluate the two-body terms 



VaiA 2 a 3 A4 — (AiA 2 |V|A 3 A4). 



(36) 



In our system the mutual interaction between the elec- 
trons is taken to be purely Coulombic, i.e. 



v(fy) 



(37) 



where e is the effective dielectric constant of the mate- 
rial. The interaction operator is thus diagonal in spin. 
Recalling that our single-particle states |A) were two- 
component spinors (jHJ, the two-body term l|36|) consists 
of sum of four terms and is of the form 



1 A2 A3 A, 



= W drd?^{r)^ * (7)V(r, ?)<tf (f) 

a. a' J 

where the summation indices take values f and J.. Fur- 
thermore, since we expressed the spatial components <\f 
as superpositions of functions g n i e [Eq. (|17|l ] , 
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E< 

n=0 



1,1 




we are ultimately led to evaluate the Coulomb matrix 
elements in the oscillator wavefunction 



e- x ' 2 x^' 2 LW{x)e m 



basis. These matrix elements can be expressed in terms 
of finite sums as 3 



•A.n 1 n 2 n 3 n 4 — w ri2^2 I i — i I M 'ii3<3 w n4i 4 / 



1+& 



V2e 2 r 



mi 



(m + |4|)! 



n 2 ! 


2 


n 3 ! 


2 


^4! 


L(«2 + N)!. 




L(«3 + N)L 




L(«4 + N)U 



m n,2 ^3 ^4 



x E E E E [«i+«4+i(i^ii + i4i-*)]! 

k;i =0 K2— k,3 — K4— 
X [K2 + «3 + 2(N + |4|-fc)]! 

x (-i) B1+B4 K + I4I)!K + |4I)! 

X Ki !k 4 ! (ni-/si)!(|^i| + Ki)!(n4-K4)!(N + K4)! 
X (-1) K2+K3 (n 2 + |^|)!(n 3 + |4l)! 

K 2 !k 3 ! (n2-K2)!(|^2| + K2)!(n 3 -/C3)!(|4|+K3)! 
[«1 + K4 + + |4| + fc)]! 

- \ Kl + k 4 + |(|4| + |4| - fc) - «]!(fc + s)! 

[«2 + «3 + |(|4| + |4| + *)]! 



K14 



s=0 

«23 

§ [«2 + ^3 + |(|4| + |4| — k) — t] l(k + t)\ 
(-l) s +* T{k + s + t+\) 



s\t\ 



2k+s+t+l ' 



(39) 



where K14 = «i + k 4 + |(|4| + |4| — &)) ^23 = «2 + 
k 3 + 3(|4| + |4| - k) and k = |4 - 4| = |4 - 4|- 
Numerical techniques to evaluate these two-body terms 
and to diagonalize the resulting Hamiltonian matrix is 
described in the Appendix. 



IV. RESULTS AND DISCUSSIONS 

For numerical evaluation of the energy spectra and the 
optical absorption spectrum for QDs with a few inter- 
acting electrons, we have considered the InAs, InSb and 
GaAs quantum dots. Parameters of these systems are 
already given in Sect. II. The energy spectra and the op- 
tical absorption spectra for these systems are described 
in the subsections below. 



A. InAs quantum dots 



FIG. 10: Same as in Fig. |HJ but for a three-electron InAs 
quantum dot. 



Our numerical results for energy spectra and absorp- 
tion spectra (dipole-allowed) for 2-4 electrons are pre- 
sented in Figs. 171121 and for various values of the SO 
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FIG. 11: Same as in[7|but for a four-electron InAs quantum 
dot. 




two-electron state | , A^ 2 ) = |Ao,A_i), where |A^) is a 
spinor with j x = t x + 1/2 = 1/2, d* 1 = 0, and |A^ 2 ) a 
spinor with j2 = —1/2 and u^ 2 — 0, i.e., both electrons 
have zero orbital angular momenta with opposite spins 
(corresponding to the spinors and 1 in Table Q] with 
J = ji +j 2 = 0) . When we increase the magnetic field the 
spin triplet configuration will become, due to the interac- 
tion, energetically more favorable. If the Lande' ^-factor 
is negative then the electrons would like to occupy states 
with orbital angular momenta and —1 with both spins 
up (i.e., states and 2 of Table[IJ). In the spinor picture 
this means that |Af 2 ) still has £2 = — 1 (J = 0) but now 
u^ 2 7^ and d^ 2 = 0. The SO interaction mixes these 
two configurations which results in a level repulsion. On 
the other hand, when the strength of the SO coupling is 
further increased, the relative significance of the Zeeman 
contribution to H. decreases. The energy shifts to states 
with J 7^ will then become energetically feasible and 
we again have crossings of levels. For increasing number 
of electrons in the dot, the energy spectra is more dense 
and exhibit additional level crossings (Figs. 151- H*2*jl . As 
a consequence, the ground state angular momentum also 
changes more frequently as compared to that of the two- 
electron case. It should be pointed out that in many- 
electron dots these level crossings and repulsion are to 
be attributed, at least partly to the mutual Coulomb in- 
teractions. The level crossings/repulsions we saw earlier 
[e.g. levels 1 and 2 in Fig. El (a)] in the single-particle 
particle picture are due to the Zeeman splitting whereas 
in interacting systems crossings occur even in the limit 
of vanishing Zeeman coupling. In the present InAs dot 
the Coulomb interaction brings, for example the singlet- 
triplet transition to much lower magnetic field (B w 2T) 
as compared to the field required in a noninteracting sys- 
tem (B w 4T). 

At moderate SO coupling strengths the absorption 



FIG. 12: Same as in Fig. [H] but for a four-electron InAs quan- 
tum dot. 



coupling strength (in meV.nm), a. As in the case of 
the non-interacting electron system, we have consid- 
ered the following parameters for the InAs quantum dot: 
m*/m = 0.042, e = 14.6, g = -14 and hu = 7.5 meV. 

A striking feature visible in the absorption spectra 
(Figs. |H1 E3 an( l EH is the appearance of discontinu- 
ities, anticrossings and new modes in addition to the two 
main (a = 0) absorption lines. These optical signatures 
of the SO interaction are consequences of the multitude 
of level crossings and level repulsions that occur in the 
energy spectra (Figs. [7| |5J and ITT|l . The latter ones 
can be attributed to an interplay between the SO and 
Zeeman couplings. In order to understand their origin, 
let us first examine the case of the two-electron system 
(Figs. l7l5|) . In our spinor notation the main contribution 
to the ground state at zero magnetic field comes from the 
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FIG. 13: Same as in but for a two-electron InSb quantum 
dot. 
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FIG. 14: Same as in Fig.[S] but for a two-electron InSb quan- 
tum dot. 



FIG. 16: Same as in Fig. |H| but for a three-electron InSb 
quantum dot. 



spectra do not essentially differ from the single-particle 
spectrum. But when the coupling strength increases the 
deviation from the pure parabolic confinement also in- 
creases which in turn implies that the lowest final states 
of dipole allowed transitions are not any more achievable 
by adding Ml ± \huj c to the initial state energies. In 
particular, this results in discontinuities and anticross- 
ing behaviors as well as appearence of new modes. As 
an illustration, let us consider the absorptions that at a 
magnetic field of B — IT take the two-electron system 
from the ground state to excited states. In the absence of 
the SO coupling the ground state is a spin-singlet state 
S = with total angular momentum J = 0. According 
to the dipole selection rules absorptions cause transitions 
to states J — ±1 and 5 = with energies AE± above 
the ground state. In Fig. [SI (d) , we note that in addition 



to the two main lines there are now two additional lines 
(at around B = 1 T) of appreciable intensity at the SO 
coupling strength a — 40 (mev.nm). Further analysis re- 
veals that the ground states still have J — and that the 
expectation value of the spin z-component is (a z ) = 0. 
The excited states also have J = ±1, as before. However, 
the final spin states can no longer be classified as singlets: 
the expectation values (a z ) vary between —0.03 and 0.39. 
When the number of electrons increases the number of 
these additional modes also increases but at the same 
time the relative intensities decrease (at each B we have 
normalized the total intensity to unity). On the other 
hand, the discontinuities as consequences of deviations 
from a parabolic confinement become more pronounced 
fFigs. I§l- I12[l . This is because there are higher angular 




FIG. 15: Same as in0but for a three-electron InSb quantum 
dot. 



FIG. 17: Same as in but for a four-electron InSb quantum 
dot. 
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FIG. 18: Same as in Fig.|S] but for a four-electron InSb quan- 
tum dot. 



FIG. 20: Same as in Fig. [8] but for a two-electron GaAs 
quantum dot. 
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Figs. 1131 1151 1171 and for various values of the SO coupling 
strength a. The corresponding optical absorption spec- 
tra are presented in Figs. El El and 1 181 As compared to 
the spectra of InAs dots a clear difference is the almost 
total absence of anticrossings and discontinuities. This 
is partly due to the very large Zeeman coupling which 
practically nullifies the SO interaction at the coupling 
strengths a we are concerned with. Another reason is the 
large kinetic energies due to the very small electron ef- 
fective mass. Because the strength of the Coulomb inter- 
action is somewhat smaller than in InAs (ei n sb > einAs) 
correlations caused by the mutual electronic interactions 
are effectively much smaller in InSb than in InAs. For 
the exploration of the SO coupling via absorption spec- 
troscopy, InSb quantum dots do not seem to be a very 
promising system. 



FIG. 19: Same as in|7|but for a two-electron GaAs quantum 
dot. 



C. GaAs quantum dots 



momenta involved in the dipole transitions. As a conse- 
quence of this the upper absorption branch now exhibits 
a rich structure while in the single-particle picture it is 
practically featureless. 



B. InSb quantum dots 

As mentioned above, in addition to the InAs quan- 
tum dots, investigation of InSb quantum dots are also 
thought to be interesting, particularly in the context of 
SO coupling effects due to the large values of \g\ and a 
[4^ | . We have considered the following parameters for the 
InSb quantum dot: m*/m = 0.014, e = 17.88,3 = -40 
and Hluq = 7.5 meV. The energy levels for InSb quantum 
dots containing 2-4 interacting electrons are plotted in 



The results for GaAs quantum dots, ones that are most 
intensely explored in the absence of SO coupling, are pre- 
sented here primarily as an academic interest. The pa- 
rameters that we have used here are: m*/mo = 0.063, e = 
12.9, g — —0.44 and Huq = 7.5 meV. Clearly, the very low 
value of the |g|-factor perhaps makes the GaAs QDs un- 
suitable for any observable effect due to the SO coupling. 
Interestingly, however, among all the three types of QDs 
studied here for optical absorptions, GaAs QDs show the 
most spectracular effects for large values of a. The energy 
levels for GaAs quantum dots containing 2-4 interacting 
electrons are plotted in Figs. El El and 1231 for various 
values of the SO coupling strength a. The corresponding 
optical absorption spectra are presented in Figs. 1201 1221 
and 1241 respectively. As mentioned above, the only ob- 
served value of a for GaAs QD reported as yet is a ~ 6 
meV.nm pfij . Reversing the arguments presented in the 
previous subsection, i.e. in GaAs a very small Zccman 
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coupling and a rather large effective electron mass but 
practically equal strength of Coulomb interaction help 
us understand why the absorption spectra of our GaAs 
dots exhibit a remarkably rich structure as opposed to 
those of the InSb and InAs dots. Finding an appropriate 
set up to generate a large a for GaAs quantum dot would 
be a major (but worthwhile) experimental endeavor. 



V. A BRIEF REVIEW OF EARLIER 
THEORETICAL WORKS 

In this section, we present a critical review of ear- 
lier theoretical reports on how the Bychkov-Rashba spin- 
orbit coupling in parabolic quantum dots were treated 

§M> M EH Ea E3, Eil Q m, HI, HE HE H3, 
I29L l30| . There are quasi-exact solutions available 
for electrons confined in a parabolic quantum dot in the 
presence of the SO interaction, but without the inter- 
electron interaction , and exact analytical results are 
also reported in the case of a circular quantum dot with 
hard walls 20], again for a non- interacting system, but 
with the Bychkov-Rashba spin-orbit interaction included. 
However, for the realistic systems of parabolic quantum 
dots with interacting electrons these methods are pro- 
hibitively complicated and evaluation of the energy spec- 
trum can only be done numerically. Among the theoret- 
ical papers dealing with the SO interaction in quantum 
dots discussed below, Kuan et al. presented the best 
treatment of the single-electron states. They looked at 
the energy levels of parabolically confined quantum dots 
with Bychkov-Rashba SO coupling and in the presence 
of zero and nonzero magnetic fields. They solved the 
single-particle equation correctly by expanding the solu- 
tion spinors in terms of the eigenfunctions of QDs with- 
out the SO interaction, i.e. the Laguerre functions. Wc 
have used a similar approach to construct the basis states 
for our multi-electron QDs [Sects. II, III]. 

Voskoboynikov et al. 0] studied the effect of SO in- 
teraction on the energy spectrum of cylindrical semicon- 
ductor QDs in an externally applied magnetic field. They 
considered the Bychkov-Rashba SO coupling due to the 
parabolic confinement, i.e. an in-plane field conserving 
orbital and spin angular momenta. As a consequence, 
the SO coupling has no qualitative effects on, for exam- 
ple, the absorption spectra. The electron-electron inter- 
action was not included in this scheme. In Ref. 16], they 
studied the magnetization and magnetic susceptibility in 
few-electron parabolic QDs with SO coupling. As in their 
earlier paper, they handled the SO term only due to the 
parabolic confinement and therefore the single-particle 
Hamiltonian is diagonal in spin space. They neglected 
the mutual electronic Coulomb interaction. 

Governale ^3] investigated the effects of SO coupling 
on the addition energy and on the spin properties of few- 
electron QDs in the absence of the external magnetic 
field. He introduced the SO coupling into single-particle 
states perturbatively but also compared the resulting en- 
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FIG. 21: Same as in0but for a three-electron GaAs quantum 
dot. 



ergies to the ones obtained by numerical diagonalization 
technique. At the small SO coupling strength that he 
considered, the perturbation approach seems to be valid. 
Electron correlations were handled by using the spin- 
density functional approach. Cremers et al. studied 
conductance and its fluctuations in the presence of SO in- 
teraction, Zeeman coupling, externally applied magnetic 
field, in a (single-electron) QD. They solved the single- 
particle equation applying an approximate unitary trans- 
formation which in leading order takes the Hamiltonian 
to a diagonal form in spin space. 

Valin-Rodriguez [2lJ considered a single electron in a 
parabolic QD with Bychkov-Rashba SO coupling. He 
performed a unitary transformation to transform the 
Hamiltonian in spin space to a diagonal form up to second 




FIG. 22: Same as in Fig. |S| but for a three-electron GaAs 
quantum dot. 
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order in SO and Zeeman coupling parameters. He showed 
that the effective SO interaction is influenced by the in- 
terplay between Zeeman and SO couplings. In Ref. p^ . 
Valin-Rodriguez et al. introduced a spatially modulated 
(in radial direction) Bychkov-Rashba coupling in single- 
electron (disk) QDs. They solved the two-component 
spinor equation and numerically evaluated the spin den- 
sity. They concluded that it is possible to confine elec- 
trons spatially with appropriate structural modulation. 
These authors also investigated the SO couplings in de- 
formed parabolic quantum dots (23|. They solved the 
single-particle equations using the approximate unitary 
transformations mentioned above. They were interested 
in the effects of spatial deformations to the spin splitting 
oscillations. They estimated the Coulomb interaction 
contribution using the time-dependent local-spin-density 
approximation [24j . 

Lucignano et al. [2f| studied the few-electron QDs in- 
cluding the mutual electron-electron interaction and un- 
der the influence of an externally applied magnetic field. 
They applied an exact diagonalization method (but with 
a rather restricted basis: 28 single-particle states deduc- 
ing from their earlier paper |43| . They particularly looked 
at the possibility to use the SO coupling to control the 
excitations under the magnetic fields which polarize the 
ground state, i.e., close to the final single-triplet transi- 
tion. The SO coupling is included in the many-electron 
Hamiltonian, but not in the basis states. They evalu- 
ated the dipole matrix elements for absorption from the 
ground state to the lowest dipole-allowed excited state. 
They claimed that there is an increase in intensity close 
to the transition to the fully polarized ground state. 

Destefani et al. |26| reported numerical results for 
energy levels and spin polarizations for one and two- 
electron parabolic QDs under a magnetic field and the 
SO coupling. For the two-electron system, Coulomb in- 
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23: Same as in0but for a four-electron GaAs quantum 




FIG. 24: Same as in Fig. |§1 but for a four-electron GaAs 
quantum dot. 



teraction is also included. They do not construct the 
correct single-electron states in the two-electron QD, the 
SO coupling is, in fact, taken into account only in the 
many-electron Hamiltonian. Debald et al. H3| studied 
oscillations in few-electron parabolic QD in a magnetic 
field, between states where the degeneracy is lifted by 
the SO coupling, i.e. at the level repulsion points. The 
Coulomb interaction was taken into account only approx- 
imately, because the many-body effects were claimed to 
play only a minor role in the very small magnetic field 
considered in that work. Konemann et al. con- 
sidered the SO coupling in single electron QDs. They 
showed that there is an anisotropy between spin split- 
tings due to magnetic fields parallel and perpendicular 
to the dot. The anisotropy was shown to be propor- 
tional to the strength of the SO coupling. Bcllucci and 
Onorato (2i| studied the influence of SO coupling on the 
charge and spin polarization in a vertical disk-shaped 
QD under a strong perpendicular magnetic field. They 
treated the SO coupling perturbatively (upto second or- 
der). They handled the Coulomb interaction within the 
Hartree-Fock approach. They studied the energy split- 
tings due to the SO coupling. Finally, Fransson et al. 
|30| studied transport through QDs with spin dependent 
couplings to the contacts. They evaluated the QD en- 
ergy levels using a (first principles) density functional 
theory. They calculated the transport properties of (a) 
non-interacting electrons taking into account the few lev- 
els closest to the Fermi level, and (b) interacting electrons 
using an approximate Hamiltonian with the levels closest 
to the Fermi level. We would like to note here that, in 
the light of all these theoretical approaches, our method 
of including the SO coupling for interacting electrons in 
a parabolic QD seems to be the most accurate one. Al- 
though, given the fact that our approach involves ex- 
tensive numerical computations, some of the approaches 
discussed above, such as the one by Lucignano et al. (2j| , 
seem to be very promising. 



16 



VI. CONCLUSIONS 

In conclusion, we have studied the energy levels and 
optical-absorption spectra for parabolic quantum dots 
containing upto four interacting electrons, in the pres- 
ence of spin-orbit coupling and under the influence of 
an externally applied, perpendicular magnetic field. We 
have presented a very accurate numerical scheme to eval- 
uate these quantities. We have presented results for the 
Fock-Darwin spectra in the presence of SO coupling for 
quantum dots made out of three different semiconduc- 
tor systems, InAs, InSb, and GaAs. The effects of SO 
coupling on the single-electron spectra are primarily to 
lift the degeneracy at B = 0, rearrangement of some of 
the energy levels at small magnetic fields, and level re- 
pulsions at high fields. These are explained as due to 
mixing of different spinor states for increasing strength 
of SO coupling. As a consequence, the corresponding 
absorption spectra reveal anticrossing structures in the 
two main lines (a — 0) of the spectra. For the interact- 
ing many-electron systems we observed the appearence 
of discontinuities, anticrossings, and new modes that ap- 
pear in conjunction with the two main absorption lines. 
These additional features arise entirely due to the SO 
coupling and are a consequence of level crossings and 
level repulsions in the energy spectra. An intricate inter- 
play between the SO coupling and the Zeeman energies 
are shown to be responsible for these new features seen in 
the energy spectra. Our accurate results for the low-lying 
energy levels for the SO coupled QDs can also be mea- 
sured, in principle, by transport [33] or capacitance (3^ 
spectroscopy, which have been successfully employed ear- 
lier to map out the energy spectra of parabolic quantum 
dots. Optical absorption spectra for all three types of 
quantum dots containing a few interacting electrons that 
are studied here show a common feature: new modes ap- 
pear, mostly near the upper main branch of the spectra 
around 2 tesla that become stronger with increasing a. 
Among the three types of systems considered here, op- 
tical signatures of the SO interaction is found to be the 
strongest in the absorption spectra of a GaAs quantum 
dot, but only at very large values of the SO coupling 
strength, and appears to be the weakest for the InSb 
quantum dots. Experimental observation of these new 
optical modes that appear solely due to the presence of 
SO coupling would be very exciting because that would 
be a major step forward in our quest to manipulate the 
spin dynamics in nanostructured systems via the SO cou- 
pling. 

Our future works along this line will be to explore cou- 
pled QDs, or QD molecules [44j. Our primary goal will be 
to generate accurate results for energy levels and optical 
absorption spectra for coupled (laterally or vertically 
|46|) quantum dots with spin-orbit interaction. In addi- 
tion to being important for fundamental studies, these 
results would be interesting from the point of view of 
quantum computations |47| as well. It is now well recog- 
nized that semiconductor quantum dots have the poten- 



tial to become the building blocks for solid state quantum 
computation. Quantum states of single, double and even 
triple coupled [48j quantum dots have been explored for 
this purpose. In a quantum computer, information is 
stored in a two-level system. Hence a promising candi- 
date system to realize the quantum bits, the fundamental 
unit of information in a quantum computer, is a quantum 
dot where the single-electron states can be used for that 
purpose. In a magnetic field, the Zeeman splitting of elec- 
tron spin can provide a two-level system (for an odd num- 
ber of electrons). Alternatively, the spatial wave function 
of a single-electron state in a double quantum dot (allow- 
ing for electron tunneling between the two dots) can also 
represent a two-level system. Spin-orbit coupling in the 
coupled quantum dot systems could perhaps be used to 
perform quantum computation (using the spin rotation, 
for example) . Accurate results for the energy levels of the 
coupled-dot system might be beneficial in that direction 
of research. The effects of SO coupling on the energy 
levels and absorption spectra for a more complex system 
such as coupled QDs are, however, important and inter- 
esting in their own ways. These would be the subjects of 
our future publications. 
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APPENDIX: DIAGONALIZATION OF MONSTER 
MATRICES 

We have mentioned in Sect. Ill about the challenging 
task of construction of the Hamiltonian matrix and its 
numerical diagonalization. Here we present a brief dis- 
cussion about the numerical method that we believe the 
most accurate (and appropriate) for that task. While 
in principle, evaluation of Eq. (|39f) is straightforward it 
turns out to be numerically highly unstable, primarily 
due to the expansion (|17f) extending to Laguerre polyno- 
mials of large degree and large angular momenta which in 
turn, leads to large terms of alternating sign |49| . A rem- 
edy for this is to employ multiple precision arithmetics 
such as, for example implemented in the Gnu arbitrary 
precision GMP library. However, if we apply multiple 
precision arithmetic directly into the sixfold summation 
the time consumed to evaluate these becomes in- 
surmountable. To circumvent this obstacle we note that 
many terms in the sums actually depend on very few 
parameters, the range of these parameters is restricted 
and the same functional forms repeat themselves. Thus 
a natural solution is to tabulate these forms and the sub- 
sums. In our Coulomb matrix element code we used the 
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tabulated functions 
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Now the summations S in the expression (|39fl can be 
written as 
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Although the summation here is still fourfold it is never- 
theless several orders of magnitude faster than the origi- 
nal one i|39|) and as such fast enough for our purposes. 

Under the influence of the SO coupling the total spin 
S of our many-electron system is not a conserved quan- 
tity. As a consequence of this we cannot fix the total S z 
of the many-body basis. This degree of freedom tends to 
make the number of non-interacting many-body states of 
the basis very large even for a small number of electrons, 
and even if the conservation of the total angular momen- 
tum J z — L + S z is taken into account. For example, to 
achieve a convergence for a four electron system in the 
parabolic confinement with harmonic potential (Hu>o) of 
few meV the size of the basis must be of the order of a 
million. Furthermore, since we want to study properties 
of the eigenstates, such as polarization and dipole matrix 
elements between states, we also need the relevant eigen- 
vectors. Clearly the sheer size of the matrix prohibits a 
full diagonalization and we have to resort to an iterative 
scheme aimed to search a given number of energetically 
lowest eigenvectors. Of course the algorithm should be 
fast and hopefully also robust. 

The algorithm proposed by Davidson and Liu (DL) 
to evaluate the eigenvalues and eigenvectors of "monster 



matrices" |5Cj seem to fit our criteria. Like any other 
iterative method it transforms the diagonalization of a 
matrix A to minimization of the Rayleigh quotient 



Ax 



\==-f±, (A.2) 
x 1 x 

where x represents the column vector of the coefficients 
in the superposition of the basis states. Also, as in many 
other methods, the only operations involving the matrix 
A are vector multiplications. This allows us to exploit 
fully the sparseness of A, i.e. we have to store only the 
non-zero elements. 

The key idea behind practically any iterative method 
subjecting the matrix only to multiplication is to search 
the minimum of the quotient (|A.2fl in a (very small) sub- 
space and to update this subspace in each iteration step. 
How this updating is performed depends on the method. 
For example, in the common conjugate gradient method 
the subspace is two-dimensional and spanned by the gra- 
dient g at current position x and a vector s conjugate to 
it with respect to A (i.e. s T Ag = 0). 

In the DL method the dimension of the search space 
varies from step to step. Suppose that at a given step 
our search space S is spanned by the orthonormal vectors 
si, . . . , sk (the dimension K must, of course, be greater 
than the number of required eigenstates). Finding the 
minimum of l|A.2|l in this subspace corresponds to the 
diagonalization of the K x K matrix S T AS (we take S 
to represent also the matrix with columns Sfe). As a result 
we get K eigenvalues /ik and K eigenvectors Zk, each of 
dimension K. The expanded vectors 



Vk 



K 

E 

i=l 



will approximate the eigenvectors and the quantities [i^ 
the eigenvalues we are seeking for. The next task in the 
iteration step is to update the subspace S. For that pur- 
pose we pick up a certain number (a parameter depending 
for example on the size of the computer memory) of the 
residuals 

r k = (A- ^ k )Vk 

with largest norms. The selected residuals are orthonor- 
malized with respect to the space S and then appended 
to it. So, in each step the dimension of the search space 
and the size required to store it increases and we may 
eventually exhaust all the memory. At this point we 
compress the space S to its bare minimum comprising 
only as many vectors as we are required to find. These 
vectors are selected from the set {yk} and are the ones 
with smallest fj,k- 
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